Localization protected quantum order 
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Closed quantum systems with quenched randomness exhibit many-body localized regimes wherein 
they do not equilibrate even though prepared with macroscopic amounts of energy above their ground 
states. We show that such localized systems can order in that individual many-body eigenstates 
can break symmetries or display topological order in the infinite volume limit. Indeed, isolated 
localized quantum systems can order even at energy densities where the corresponding thermally 
equilibrated system is disordered, i.e.: localization protects order. In addition, localized systems 
can move between ordered and disordered localized phases via non-thermodynamic transitions in 
the properties of the many-body eigenstates. We give evidence that such transitions may proceed 
via localized critical points. We note that localization provides protection against decoherence that 
may allow experimental manipulation of macroscopic quantum states. We also identify a 'spectral 
transition' involving a sharp change in the spectral statistics of the many-body Hamiltonian. 



Our current understanding of the phases of quantum 
matter in equilibrium is built largely on the traditional 
Landau framework of broken symmetries [1 and the more 
recent framework, still in rapid evolution, of topological 
order and allied classifications [2]-[7]. There are inter- 
esting exceptions to these in the presence of quenched 
randomness: such as Anderson localization [8 , which 
is firmly established in studies of non-interacting parti- 
cles [9^. Recently, the work of Basko, Aleiner and Alt- 
shuler [10 and others pTHTS] has added to these a pre- 
viously conjectured [8] extension of Anderson localiza- 
tion to closed, quantum interacting systems — the phe- 
nomenon now known as many-body localization (MBL). 

To understand the nature of many-body localization, 
it is useful to first refer to the eigenstate thermalization 
hypothesis (ETH) p^HlS]. The ETH, when true, ap- 
plies to the exact many-body eigenstates of the Hamil- 
tonian of a closed, isolated quantum system, in the limit 
of many degrees of freedom. The ETH postulates that 
for a large class of quantum systems, the probability op- 
erator (a.k.a. the reduced density matrix) for any sub- 
system is, in any exact many-body eigenstate of the full 
system, equal to the equilibrium Boltzmann-Gibbs dis- 
tribution at the temperature set by the energy density 
of the eigenstate. This occurs because the remainder of 
the full system successfully acts as a thermal bath for the 
subsystem in question. 

The many-body localization phase transition is an 
eigenstate phase transition from the thermal phase where 
the exact many-body eigenstates obey the ETH, to the 
localized phase where the eigenstates violate the ETH; 
the latter fail to be a heat bath that thermally equi- 
librates its subsystems [HI [TOl IHl US]. Thus this is a 
dynamic, but not thermodynamic, phase transition from 
the thermal phase where the system does thermally equi- 
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FIG. 1: In the one-dimensional models we consider, eigen- 
state properties depend on the energy density, the typical 
value of /i/J, the disorder strengths 5j and 5h^ and the 
fermion-fermion interaction strength A. This figure shows, 
schematically, a slice through the phase diagram, for typical 
h/ J < 1, and non-zero A. Extended (Ext) and many body lo- 
calized (MBL) phases are separated by the localization transi- 
tion, ordered spin glass (S-G) and paramagnetic (Para) phases 
are separated by an eigenstate phase transition, and regions 
with paired and unpaired spectra are separated by a spectral 
transition. 



librate under the dynamics due to its own Hamiltonian, 
to the 'glassy' localized phase where the isolated quan- 
tum system can remain far from thermal equilibrium for- 
ever. In a related diagnostic, MBL eigenstates generally 
display an 'area law' entanglement entropy [19 , unlike 
thermal eigenstates where the entanglement entropy is 
generally a 'volume law' that reproduces the thermody- 
namic entropy. 

In this work we examine the highly excited eigenstates, 
defined as eigenstates with a macroscopic energy above 
the ground state, of MBL systems and point out that 



they come in many flavors, and may be classified in terms 
of broken symmetries, topological order and/or critical- 
ity, very much as in the usual account of phases and 
phase transitions in equilibrium systems. We note that 
in the presence of many-body localization, equilibrium 
constraints on order can be evaded: symmetry breaking 
can occur in highly-excited states of one-dimensional sys- 
tems, and topological order can arise even in the absence 
of a hulk gap. Instead it is the localization that 'pro- 
tects' the topological order. We suggest that there can 
even arise continuous phase transitions between distinct 
MBL phases, which proceed via MBL critical points. Fi- 
nally, we point out the existence of yet another kind of 
transition within the MBL phase: a 'spectral transition', 
which does not involve a change in the properties of the 
eigenstates, but instead a change in the spectral statistics 
of the system's Hamiltonian. We emphasize that unlike 
standard discussions of quantum phase transitions, our 
discussion is not about ground states or low-lying excited 
states, but is about highly-excited eigenstates at energies 
that would correspond to nonzero (even infinite) temper- 
ature if the system could thermalize at these energies. 
This might be useful for experiments, in particular for 
experiments designed to exploit topological order. 

This article is structured as follows: We first con- 
sider three equivalent "integrable" [20] one-dimensional 
systems - the Majorana chain, p-wave superconducting 
chain, and transverse-field Ising chain [21 . The eigen- 
states of the disordered versions of these Hamiltonians 
are localized and can be described in terms of noninter- 
acting localized fermions. This localization is robust to 
weak fermion-fermion interactions ^ . The MBL eigen- 
states can break symmetries, and a robust notion of topo- 
logical order can be defined for them even without a gap. 
We also show how by tuning parameters of the Hamil- 
tonian, we can drive an eigenstate phase transition from 
one MBL phase to another. In the Ising chain, the transi- 
tion is of the symmetry-breaking type, from a disordered 
paramagnet to a phase where the eigenstates are feline 
(Schrodinger cat) states with long-range spin-glass (SG) 
order. In the Majorana and Dirac fermion systems, the 
transition is between MBL states with and without topo- 
logical order. We further argue that the critical point 
separating these two distinct MBL phases can itself be 
MBL. We close by discussing extensions of these ideas to 
higher dimensions. We emphasize again that the ideas 
discussed here apply only to closed quantum systems - 
i.e. it is essential that the system not be put in contact 
with an external thermal bath. 



The models: The three non-interacting one- 
dimensional models studied in this paper (to which 
we will add interactions perturbatively), are the trans- 
verse Ising, Majorana, and Dirac fermion chains of finite 



length L sites, with open ends: 
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where the a^'^ are Pauli matrices, the a, h are (self- 
adjoint) Majorana fermion operators and the c, c^ are 
conventional Dirac fermion operators. The parameters 
Ji and hi are drawn from distributions P{J) and P{h) 
with means J > and /i > 0, and variances ^j and 5\. 
The precise details of the distribution are unimportant 
for our present purposes, but it is vital that at least one 
of the variances be non-zero. For specificity, take the dis- 
tributions to be log-normal, so all Ji and hi are positive. 
We consider the disorder strengths 5'^ to be parameters 
in our analysis. These three Hamiltonians are related to 
one another by the duality transformations [2TJ |22] 
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The above Hamiltonians all possess a global Z2 symme- 
try, implemented by the operator P, which takes the form 
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In the ordered phase (typical h/J < 1), there is an 
'edge mode' operator O^, whose commutator with the 
Hamiltonian is exponentially small {{h/J)^ <C 1) in L. 
In the Majorana language [21], 
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This operator creates a Dirac fermion which is bilocalized 
around the two edges, and has an energy which is expo- 
nentially small in the system size, E ^ exp ( — L ln( J//z)) . 
The existence of this quasi-zero-energy edge mode can 
be considered a diagnostic for topological order. Since 
O and P anti-commute, acting on a state with O or O^ 
flips its eigenvalue under P. 



Localization in the ordered phase: To discuss lo- 
calization in the non-interacting model, it is convenient 
to use the Ising spin formulation ([l]). The results apply 
equally to the Majorana and Dirac fermion chains. 

In the ordered phase, J — V^ ^ /i, the ground state 
is a ferromagnet, with long-range order of a^. It is even 
under P, consisting of a linear combination of states with 
the average z-magnetization 'up' and 'down'. Let's call 
it |0, +). It is a 'feline' (Schrodinger cat) state consisting 
of a coherent linear combination of two 'macroscopically' 
different states. If we let O operate on this ground state, 
this makes the other 'ground state', |0, — ) = O|0, +), 
which is also feline, is odd under P, and is higher in 
energy by an amount that is exponentially small in L. 

The normal modes in the ordered phase are fermion 
operators, which anticommute with each other and with 
P. The action of these normal modes on the spin state is 
straightforward: they create domain walls (with respect 
to the ferromagnetic ground state). Domain walls sit on 
bonds, and are created by the self-adjoint operators 
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(7) 
In the integrable Hamiltonian, the properties of an 
eigenstate may be extracted from the behavior of the 
normal modes, which create domain walls. In the clas- 
sical limit /i = 0, the normal modes are trivially local- 
ized. Away from the classical limit, they hop under the 
action of the transverse fields /i^, and see a spatial po- 
tential Ji. In the clean limit, the normal modes (and 
hence the domain walls) are delocalized over the chain. 
However, it is well known that non-interacting fermions 
hopping in a one-dimensional random potential always 
experience localization [9]. Therefore, in the presence 
of non- vanishing disorder in J^, the normal modes (and 
hence the domain walls) must necessarily be localized. 

2 

In the strong disorder limit, when (5j ^ /i , the single 
domain wall eigenstates are exponentially localized, with 
localization length 
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This follows straightforwardly from perturbation theory 
in small h/3j. Localization of the domain walls in each 
many-body eigenstate implies that each eigenstate has 
long-range spin glass order [23 , i.e. the correlation func- 
tion {(jQCFr) within the eigenstate is non-zero for large 
r, with its sign set by how many domain walls are 
present and localized between sites and r. Observe 
that this spin-glass order breaks the global Z2 symme- 
try in these highly-excited localized eigenstates of this 
one-dimensional system, although at thermal equilibrium 
such discrete symmetry breaking at non-zero tempera- 



tures is forbidden. This eigenstate spin-glass ordering is 
also discussed in Pekker, et al. j24j. 

Just as this model in this ordered phase has two nearly- 
degenerate ground states |0, ±) that are each feline, the 
excited eigenstates also come in nearly-degenerate fe- 
line pairs |n, ±), produced by adding some particular 
set of localized domain walls, n, to each of the two 
ground states. The states |n, ±) are orthogonal eigen- 
states of the global spin fiip operator P, with eigenvalue 
±1. From these two eigenstates we can make the state 
|n,t) = (1^7+) + 1^5 ~))/v^ which has some particular 
pattern of local z-magnetizations dictated by the local- 
ized domain walls. Each eigenstate |n, ±) is a (feline) 
coherent linear combination of this spin-glass state |n,t) 
with the opposite state under the global spin fiip opera- 
tion, \n^\) = P|n,t)- The energy splitting between |n, ±) 
is exponentially small in system size, and thus the energy 
uncertainty in the symmetry-breaking states |n,t) and 
|n,|) vanishes exponentially with the number of spins in 
the infinite volume limit. We note in passing that the ex- 
istence of spin glass order in the eigenstates is consistent 
with exponential localization of the local fiuctuations of 
conserved quantities - a hallmark of the MBL phase. 
Robustness of localization to weak interactions: 
We assume we are working with highly-excited states, 
where the localized domain walls are dense and frequently 
overlap (the typical domain wall separation can be less 
than a localization length) [25 . We then add to the 
Hamiltonian weak short-range interactions between these 
domain walls, and ask whether they cause a breakdown 
of localization. This situation is essentially identical to 
that analyzed by Basko, et al. [IO4, who show that many- 
body localization is stable to weak nonzero interaction A 
as long as the localization length of the single-particle 
states is finite. We also assume that the interactions 
commute with P and thus respect the Z2 symmetry. 

If the disorder is weak, the localization will be de- 
stroyed by strong enough interaction, making the ex- 
tended, thermal phase at weak disorder and high en- 
ergy density, as shown in Fig. 1. In terms of the Ising 
model, in the extended phase the spin-spin correlations 
within a thermal eigenstate are short-range. If we assume 
the eigenstate phase transitions are continuous (not first- 
order), then there must be two transitions, as in Fig.l. 
Starting in the extended thermal phase and increasing 
the disorder, first we reach the localization transition. 
Once we enter the localized phase, the spin-spin corre- 
lations in the eigenstates can start increasing to longer 
range than they are at thermal equilibrium. The transi- 
tion to the MBL spin-glass phase is when the spin corre- 
lations develop long-range order. 

Spectral transition in the ordered phase: In the 
MBL spin-glass phase, all eigenstates come in parity- 
related pairs |n, ±) which differ only by their occupation 
of the single particle edge mode created by O^. The en- 
ergy of the edge mode is exponentially small in the chain 



length, Eo ~ exp(— L/^), where ^ is the locahzation 
length of the edge mode. Thus, the level splitting within 
such a pair is Eq ^ exp(— L/^). Meanwhile, the typi- 
cal many-body level spacing is of order SE ~ exp(— sL), 
where s is the thermodynamic entropy per site that would 
result if the system equilibrated at the given energy den- 
sity. Thus, tuning ^ or s leads to a spectral transition. 
For strong disorder and thus small localization length ^ 
and/or for low energy density and thus low entropy s, the 
spectrum is 'paired'. Here the edge mode level splitting 
is less than the typical level spacing so the many-body 
spectrum at large L consists of nearly-degenerate dou- 
blets with Poisson inter-doublet level spacing statistics. 
The level spacing within each doublet is exponentially 
small in system size compared to the typical level spac- 
ing between doublets. Meanwhile, for weaker disorder 
and/or higher entropy, the spectrum is unpaired. In this 
unpaired regime, the energy of the single particle edge 
mode is exponentially larger than the typical many body 
level spacing, and as a result we have Poisson statistics for 
all of the individual many-body energy levels. The Pois- 
son (or paired Poisson) level statistics themselves are a 
diagnostic for MBL, since in the thermal phase, the level 
statistics are instead characteristic of the Gaussian or- 
thogonal ensemble (GOE). 

The spectral transition is illustrated in Fig. [2] It does 
not involve a change in the symmetry or topological or- 
der of the eigenstates, but rather involves a change in 
the many-body spectral statistics of the Hamiltonian. It 
may be detected in numerics by examining the ratio of 
two consecutive gaps, r = min(J^, ^^+i)/max((5n, (^n+i)^ 
where Sn is the energy gap between the n^^ and (n + 1)^^ 
many body eigenstates. This ratio will have average value 
(2 In 2 — 1) in the 'unpaired' localized regime [T2^, whereas 
the average value in the paired regime will be exponen- 
tially small in system size (zero in the infinite volume 
limit). Thus the ratio of two consecutive gaps provides a 
sharp diagnostic for the spectral transition. 
Magnetic response: One curious feature of these pairs 
|n, ±) of feline eigenstates is that the states |n,t) have 
a typical z-magnetization m that scales as ~ a/L. A 
longitudinal magnetic field B added to the Hamiltonian 
thus acts on the Hilbert space of a pair (in the B = 
eigenstate basis) as 
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It is easy to verify that the energy eigenstates have an 
adiabatic magnetic susceptibility x = {dm/dB)^^^ that 
is exponentially large in L, but is of opposite sign between 
the two paired eigenstates. The feline nature of the eigen- 
states is thus destroyed by an exponentially small (in L) 
longitudinal field. Interestingly, this makes this system 
a possibly very sensitive magnetic field sensor. If a state 
with a nonzero z-magnetization is prepared, it will be lo- 
calized and stable as long as the magnitude of the longi- 



FIG. 2: An illustration of the spectral transition between 
paired and unpaired many-body spectra, driven for example 
by tuning disorder strength. 



tudinal field B is larger than ~ e~^l^ . However, once the 
magnitude of the longitudinal field drops below this level, 
then the system's many-body eigenstates, although still 
localized, are feline and the magnetization will decrease. 
The main caveat to this apparently exponentially-fine- 
in-L field sensitivity is that the resulting magnetization 
dynamics is exponentially slow in L, and the utility of 
this idea will be limited by the decoherence rate of any 
real system. We note however that, in our three dimen- 
sional world, three dimensional MBL systems where only 
the boundary is coupled to the environment should be 
protected from environmental decoherence, since the ef- 
fects of the environment should only propagate into the 
system up to one localization length. 
Localization in the paramagnetic phase: We now 
turn our attention to the disordered regime, h — \fo\ ^ 
J. In this regime, the (Ising) ground state has all spins 
aligned on average along the x axis, and the elementary 
excitations are spin fiips. A similar analysis to above 
leads one to conclude that all spin fiips are localized in the 
non-interacting model, as long as J^ ^ 0. Localization is 
robust against weak enough interactions so long as the 
localization length ^ is finite. Thus, the eigenstates of the 
disordered Ising paramagnet can also be MBL. However, 
these eigenstates do not break Z2 symmetry and do not 
come in pairs. Likewise, there is no edge mode in the 
fermion language, and no topological order. 

These results can also be understood in terms of the 
well known self-duality of the one dimensional Ising 
model [22 , which swaps h and J. In particular, it follows 
from self duality that if the ordered phase is MBL, the 
disordered phase is also MBL. However, the spectrum in 
the paramagnetic phase is not paired. This is because the 
pair of parity related MBL eigenstates |n,±) on the or- 
dered side differ by the edge mode O, and the self duality 
acts non-locally on the edge mode, mapping the 'parity 
related pair' of ordered states to disordered states with 
different boundary conditions. Thus, the disordered side 



(with specified boundary conditions) has states that do 
not come in pairs and do not have edge modes, and are 
separated from their counterparts on the ordered side by 
a phase transition. 

The phase transition: In the absence of interactions 
(A = 0) there is a continuous transition between local- 
ized phases, which proceeds via an infinite randomness 
critical point even for highly excited states. We begin by 
discussing this non-interacting picture, before examining 
the effect of interactions. The critical regime is treated 
using the strong disorder renormalization group (SDRG) 
[26] . The SDRG proceeds by sequentially identifying the 
strongest bond or field in the Hamiltonian, diagonalizing 
it, and determining the coupling to the rest of the sys- 
tem perturbatively. This procedure can be shown to be 
asymptotically exact, because the system fiows to strong 
disorder. 

When looking for ground states [26], after diagonal- 
izing a particular bond or field, we should truncate to 
the lower energy subspace. However, this method can be 
easily generalized to obtain excited states, by sometimes 
truncating into the higher energy subspace instead. It 
can be verified that this does not change the fiow equa- 
tions for \Ji\ or \hi\^ but merely introduces some extra 
minus signs (some bonds become antiferromagnetic and 
some fields point along —x). We thus recover the flow 
equations [26] for the probability distributions P{h) and 
P{J) as flow equations for P{\h\) and P(| J|). The flow 
[26] is to strong disorder, so the transition should proceed 
via an inflnite randomness critical point, for all excited 
states. 

In either phase we have already argued that the non- 
interacting system is localized. However, at criticality 
{e = 0), the spectrum of the non-interacting Hamilto- 
nian contains states with all localization lengths and the 
single-fermion states in the limit of zero energy have a 
diverging localization length. Nonetheless, the critical 
point of the non-interacting Hamiltonian ( 1|2|3 ) is lo- 
calized in the following sense. Almost all single-fermion 
wave functions are localized, with only the limiting zero- 
energy states having inflnite localization length. The en- 
tanglement entropy (within a many-body eigenstate) of 
a subregion is thus sub-extensive, unlike a thermal state, 
which has extensive entanglement entropy pjj. In fact, 
since the KG flow for P(| J|), P(|/i|) is the same in the 
excited eigenstates as for the ground state, the entangle- 
ment entropy is also the same in the excited eigenstates 
as in the ground state. As determined in [27], the entan- 
glement entropy of a subregion at the inflnite random- 
ness critical point has a leading term that scales as InL, 
whereas a thermal state would have an entanglement en- 
tropy that scales as L. Thus, the (non-interacting) criti- 
cal point violates the ETH, and can be considered local- 
ized. 

We now discuss whether localization can survive inter- 
actions at the critical point. The key question is whether 



the presence of a subextensive number of critical modes 
can cause the rest of the degrees of freedom to delo- 
calize, by mediating resonances between distant near- 
degenerate localized modes. In [28 , it was determined 
that the mediated interactions fall off exponentially with 
distance, whereas the typical level spacing decreases as 
a power law of distance. Thus, the mediated interac- 
tion between distant near-degenerate modes should be 
weaker than the level splitting, and mediated interac- 
tions should be unable to delocalize the formerly local- 
ized degrees of freedom. This argument suggests that 
the critical point separating two MBL phases can itself 
be MBL. A slice through the phase diagram is presented 
in Fig[T] For more detail on the phase transition within 
the MBL regime, see [23]. 

Topological order in the fermion language: Trans- 
lating to the Majorana and Dirac fermion languages, we 
again conclude that all bulk modes are localized, with 
localization lengths equal to those calculated above. The 
bulk normal modes are Dirac fermions. Meanwhile, the 
nearly-degenerate pairs of excited eigenstates are states 
where the single particle edge mode created by O^ is ei- 
ther occupied or unoccupied. Thus, the eigenstates of 
( 2p ) have topological order. The existence of topologi- 
cal order follows trivially in these non-interacting models, 
because the edge mode does not interact with anything in 
the bulk, and is bilocalized as two Majorana modes, one 
near each end of the chain. In the clean system, the ad- 
dition of arbitrarily weak interactions destroys the topo- 
logical order, since the Majorana end modes can couple 
through the delocalized domain walls. However, in the 
disordered system, we will show that many body local- 
ization protects the topological order. 
Localization protects topological order: We now 
demonstrate that the edge Majoranas remain localized in 
the presence of interactions, even in the absence of a gap. 
We assume that the interaction locally conserves fermion 
number modulo 2, and is thus a product of an even num- 
ber of Majorana operators. Before turning on interac- 
tions, there are an odd number of quasi-zero-energy Ma- 
jorana modes localized at each edge. We cannot turn this 
into an even number of Majoranas by acting with an even 
number of Majorana operators, so the edge Majorana 
must survive the addition of interactions. The Majorana 
cannot disappear because of hybridization with localized 
bulk modes because the bulk modes are Dirac. The only 
way to make the topological order disappear is to couple 
the two Majoranas at either edge. However, the inter- 
action cannot do this because it is short range, and the 
Majorana cannot be passed from one bulk Dirac mode to 
another, until it reaches the other Majorana at the oppo- 
site end of the chain, because the hulk is MBL^ and does 
not allow energy or particles to propagate. Thus, MBL 
in the bulk protects topological order in highly-excited 
localized states, just as a bulk gap can protect topolog- 
ical order in ground states. This is in sharp contrast to 



the clean (non-MBL) system, where the Majorana edge 
modes can hybridize with each other through delocahzed 
bulk modes in any excited state. The localization protec- 
tion of the edge Major anas in quantum states other than 
ground states might be useful for experiments designed to 
exploit topological order. In particular, it might be use- 
ful for experiments designed to detect Majorana fermions 
in quantum wires [29 . 

Symmetry Breaking in d > 2: The ideas discussed 
above have a straightforward extension to the Ising model 
in more than one dimension. A major difference is that 
in two or more dimensions, thermodynamic Z2 symmetry 
breaking persists to nonzero excitation-energy densities 
even in the extended (thermal) phase. As a result, the 
phase diagram contains another type of phase transition - 
the usual thermodynamic phase transition between states 
with and without ferromagnetic order. A slice through 
the {d > 1) -dimensional phase diagram is presented in 
Fig.i 

As in Fig. 1, we assume that the eigenstate phase tran- 
sitions are continuous (not first-order). This again im- 
plies that a localized paramagnetic phase exists between 
the extended (thermal) paramagnetic phase and the lo- 
calized spin glass in the higher-energy regime above the 
ferromagnetic phases. Moving across this phase by in- 
creasing disorder, the spin-spin correlation length grows 
continuously from the finite thermal value of the ex- 
tended phase and diverges at the transition to the spin 
glass. We note that this symmetry-breaking paramagnet- 
to-spin-glass phase transition in the localized Ising model 
in d > 2 is also governed by an infinite randomness fixed 
point [30l |3T] which again should extend to finite en- 
ergy density states and lead to sub-thermal entangle- 
ment. We believe that localization persists here too, 
although a more detailed analysis is desirable. Assum- 
ing the symmetry-breaking transition in the localized 
phase out of the paramagnet is indeed in the infinite- 
randomness universality class, there should be localized 
domain walls in the ordered phase just across the tran- 
sition. The higher-energy eigenstates with these domain 
walls present are spin-glass states. This is why we believe 
the spin-glass phase always exists between the param- 
agnet and the ferromagnet in the localized regime, as 
shown in Fig. 3. We also note that unlike in (i = 1, the 
nearest-neighbor transverse-field Ising model in d > 2 is 
an interacting system on its own. 

Topological Order in d > 2: From our results thus far 
we can immediately draw some interesting conclusions 
about topologically ordered systems in d > 2 by simply 
dualizing the Ising model. This leads to MBL protected 
topological order in higher dimensions, as we now discuss. 

Let us begin in d = 2 where the "Ising model" of a 
system that exhibits topological order is the the Z2 lattice 
gauge theory with matter, governed by the Hamiltonian 
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FIG. 3: A (schematic) slice through the {d > 1) -dimensional 
phase diagram at non-zero A and typical h/J<l. In addi- 
tion to the eigenstate phase transitions and spectral transi- 
tion present in one dimension (Fig. [l]), now there is also the 
usual thermodynamic phase transition between phases with 
and without ferromagnetic order, i.e. net spontaneous mag- 
netization. The presence of ferromagnetic order is labelled by 
F. 
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supplemented by the constraint that we restrict its action 
to "gauge invariant" states defined by 
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(11) 



In the above the gauge {a\) and matter (rj) operators 
act in spin- 1/2 Hilbert spaces that live on the links / 
and sites s of the square lattice with plaquettes p and 
dp and dl are the boundaries of the corresponding ob- 
jects. On dualizing the Ising model, we get the Z2 gauge 



theory without matter, obtained from (10) by dropping 



the matter degrees of freedom entirely. In this case the 
parameters Kp and Vi are, numerically, the on-site fields 
and bond interactions of the dual Ising model. 

Let us briefly review some salient facts about the non- 
random system. In this, the paramagnetic and ferro- 
magnetic phases dualize, respectively, to the (topologi- 
cally ordered) deconfined and (non-topological) confined 
phases of the gauge theory where the terminology refers 
to the energy needed to separate two test charges to in- 
finity. For our purposes, it is more useful to consider the 
standard equal-time diagnostic which can be evaluated in 
individual eigenstates, the Wilson loop [22 for a contour 
C 



W[C] = {\[cjt) 



(12) 
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charged matter, as in (10), which is known to leave 



FIG. 4: The contours used to define the Fredenhagen-Marcu 
order parameter for a translationally non-invariant system. 



which exhibits a perimeter /area (P/A) law decay, 
log VK[C] ex — P/ — A, in the topological/non-topological 
phase. A useful picture of the ground states and excita- 
tions is obtained by thinking in the basis of eigenstates 
of Gi. At r ^ i^, the non-topological ground state has 
(Ji = 1 on all bonds. The elementary excitations are 



small loops of bonds where a^ 



-1, which we refer 



to as bonds with Z2 "electric" flux. These loops are 
dual to domain walls above the ferromagnetic state in 
the Ising language. As we pass to the topological phase 
at r <C i^, the loops/domain walls proliferate and their 
condensation signals the transition. The elementary ex- 
citations in the topological phase are visons, plaquettes 
where Hp ^i = ~1- 

Now the translation is straightforward and we conclude 
that with sufficient randomness in the couplings there ex- 
ist both MBL localized topological and non-topological 
phases in the Z2 gauge theory without matter. [In Figure 
3, we can simply relabel the paramagnetic and spin glass 
phases as topological and non-topological respectively.] 
Of maximum interest is the MBL topological phase whose 
topological order is protected by localization and would 
not exist in its absence at nonzero energy densities (the 
dual of the thermal paramagnet phase in Fig. 3 is a phase 
where topological order is destroyed by thermal fluctua- 
tions). In the MBL topological phase, which is usefully 
described as a state with a finite density of localized vi- 
sons, the eigenstates display a "spin glass" version of the 
perimeter law in which the magnitude of W[C] decays ex- 
ponentially with perimeter but with a sign that depends 
on how many localized visons are encircled by the loop. 
By contrast, in the non-topological phase, W[C] exhibits 
an area law. 

This account of the Z2 topological phase can be ex- 
tended in two directions. First, one can include gauge 



the clean system topologically ordered at T = when 
J <C r^. In the presence of sufficient randomness, the 
matter excitations will be localized and we will obtain 
an MBL protected topological phase in the presence of 
dynamical matter at nonzero energy densities. However 
this phase can no longer be diagnosed by examining the 
Wilson loop, which exhibits a perimeter law in all phases 
in the presence of matter. Instead we turn to a version 
of the Fredenhagen-Marcu order parameter [33 , which 
measures the "line tension" of Z2 electric fields specified 
by their distribution in an eigenstate. The version needed 
for our random setting is 
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where the notation is explained in Figure |4] and the ex- 
pectation values are taken in a particular eigenstate of 
the Hamiltonian. Extending the discussion in [33 to the 
present case we propose that R{L -^ 00) = in the topo- 
logical phase, but not in the confined phase. 

The second extension we propose involves two other 
common diagnostics of topological order in the clean 
system — a ground state topological degeneracy of 4^ on 
closed manifolds of genus g and a topological entangle- 
ment entropy of log 2 [34] . In the MBL Z2 topologically 
ordered phase we expect that the dominant finite size 
effect in large systems of linear dimension L will arise 
from 0{e~^) tunneling between clusters of 4^ finite en- 
ergy density states which (roughly) exhibit the same pat- 
tern of localized excitations but differ in the presence or 
absence of visons threading the non-contractible loops on 
the manifold. Likewise individual eigenstates will exhibit 
an area law in entanglement and a topological constant 
piece of size log 2 which can be detected by means of the 
subtraction procedure outlined, e.g. in |35l f36]. 

Dualizing the Ising model in d > 2 will yield a d-form 
gauge theory [40] with a topological phase with pointlike 
excitations and MBL stabilized topological order. Get- 
ting to a more conventional gauge theory in d > 2 how- 
ever requires some fresh thinking. For example, in d = 3 
we need to consider the localization of stringlike excita- 
tions (vison loops) and one cannot simply appeal to the 
body of results to date in this paper or elsewhere. 

Finally, we would draw the reader's attention to 
[SHI |39] where the dynamical localization of excitations 
in Z2 gauge theory with matter has been discussed in the 
language of the perturbed toric code relevant to quantum 
information storage. 

Conclusions and outlook: Thus we have shown that 
eigenstates of MBL systems come in many flavors, and 
may be classified in terms of broken symmetries, topolog- 
ical order, and criticality, just like extended states. Lo- 
calization itself can protect order through the intuitive 
mechanism of localizing excitations that would disrupt 



it. We have also identified a new kind of transition - 
the spectral transition, involving a change in the spec- 
tral statistics of the Hamiltonian. The protection of or- 
der and quantum coherence by localization might open 
the door to a new generation of quantum devices that are 
immune to environmental noise, and are not restricted to 
ground states or low-energy states. It may also be useful 
for ongoing experiments attempting to observe Major ana 
fermions in quantum wires [29] . 

We have focused, for pedagogical clarity, on broken 
Z2 global symmetry and Z2 topological order. Gener- 
alizations should be straightforward to other problems 
where the elementary excitations subversive of ordering 
can be localized at all energies by sufficiently random 
couplings. Immediate examples are p > 3 Zp clock mod- 
els in d > 1 and dual parafermionic systems in d = 1 
[2T] and Z^ gauge theories in d > 2. Farther afield we 
should flag other models with broken discrete symme- 
tries and also topologically ordered systems with discrete 
gapped excitations such the Levin- Wen models that ex- 
hibit non-abelian phases [37]. However, it is essential for 
our purposes that the system should not support long- 
range interactions, including those which might be me- 
diated by Goldstone modes or gapless gauge excitations 
("photons") that do not localize [41]. Thus, an extension 
of these ideas to systems with broken continuous sym- 
metries and continuous gauge groups looks problematic. 
However, an extension to continuous symmetries might 
be possible if the Goldstone mode were gapped out by the 
Anderson-Higgs mechanism, or by placing the system on 
a Bethe lattice, where Goldstone bosons are absent [42 . 
We defer further consideration of these issues to future 
work. 

While this work has focused on MBL eigenstates, an 
experimental construction of a system with a Hamilto- 
nian that displays many body localization will necessar- 
ily start with an initial state that is not an exact eigen- 
state of the Hamiltonian. We expect, based on analogies 
to the 'diagonal ensemble' viewpoint advocated in [43] , 
that for many initial states the density matrix at long 
times can be treated as being diagonal in the basis of en- 
ergy eigenstates. However, a detailed understanding of 
the dynamical evolution of an initial superposition state 
(or mixed state) is an important topic for future work. 

A final set of interesting open questions involves 
whether continuous phase transitions between distinct 
MBL phases necessarily proceed via a localized critical 
point. We have argued that for MBL phases in the one- 
dimensional random Ising model, phase transitions be- 
tween MBL states proceed via an infinite-disorder criti- 
cal point that is itself MBL. It remains to be determined 
to what extent this holds true for more complex models. 
Another open question involves the nature of the phase 
transition between extended (thermal) and MBL phases, 
which to our knowledge has not yet been determined in 
any system. 
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